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LIST  OP  SYMBOLS 


a 

One  half  the  distance  between  the  3ource  and 
sink  of  the  Rankine  Ovoid 

b 

Maximum  radius  of  body 

CK 

M 

Moaent  coefficient  »  — - - 

±f>c2rb2l 

c 

Constant  unifora  streao  velocity 

d 

Kaximun  dlaaeter  of  body 

F 

Froude  nuaber  based  on  length  = 

F 

Force  vector 

f 

Distance  of  source  or  sink  below  thu  undis¬ 
turbed  fluid  surface 

S 

Acceleration  of  gravity 

Ko 

, 

Modified  Bessel  functions  of  the  second  kind 
(Reference  13,  p.  78) 

ko 

g/c- 

*LI« 

i 

A  function  defined  on  page  15 

£ 

Over-all  length  of  body 

K. 

•V 

Moaent s 

a 

The  strengtn  of  a  source  (a  source  of  strength  a 
eaits  a  voiuae  per  unit  tise) 

p 

Placed  before  an  integral  sign  means  that  the 
Caucny  principal  value  of  the  integral  is  to 
be  taken 

Resultant  fluid  velocity  vector  at  the  location 
of  «  source  due  to  all  other  sources 

V 

1  T  ’  S: 

Magnitude  of  components  of  q  in  x,  y,  z  direc¬ 
tions.  respectively 

r 

Radial  distance 

u. 

v,  w 

Magnitude  of  components  of  local  velocity  in  x, 
y,  z  directions,  respectively 

X, 

y.  * 

Rectangular  coordinates 

a 

depth  of  centerline  of  body  _  f 
diameter  of  body  d 

P 

length  of  body  _  £ 
dlaaeter  of  body  ? 

t 

distance  between  source  and  sink  2a 
length  of  body  ~  T 
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x 

i 

ti 

p 

OJ 


A  parameter 

Christoffel  number  used  in  numerical  quadrature 
formula 

Strength  of  doublet 

Hass  density  of  fluid 
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Hartley  L.  Pond 


ABSTRACT 

T*e  noeaent  acting  on  a  Raaidne  Ovoid  moving  under  the  free  surface  ol  a 
fluid  of  infinite  depth  is  calculator.  In  order  to  account  fur  the  affect  at  the  waves 
formed  by  the  motion  at  the  body  a  correction  is  given  in  the  fora  of  a  second 
apprcximaUos. 


INTRODUCTION 

In  the  calculation  of  the  resistance  of  bodies  moving  below  a  free 
surface,  it  is  usually  assumed  that  the  singularity  system  (sources,  sinks, 
doublets,  etc.}  used  to  represent  the  body  can  be  taken,  as  a  first  approxi¬ 
mation.  to  be  the  same  as  the  singularity  system  representing  the  body  in  an 
unbounded  fluid.  In  the  present  report.  It  is  shown  that  in  calculating  the 
moment  acting  on  a  body  moving  below  a  free  surface.  It  is  necessary  to  mod¬ 
ify  the  singularity  system  to  account  for  the  waves  formed  by  the  motion  of 
the  body.  To  the  order  of  approximation  considered  in  the  present  report  the 
modification  of  the  singularity  system  is  such  that  it  gives  rise  only  to  a 
couple.  That  is,  there  is  no  change  lr.  the  horizontal  and  vertical  forces 
acting  on  the  body. 


THE  FIRST  APPROXIMATION 

With  the  usual  assumptions  that  the  wave  slope  is  small,  and  that 
the  velocity  (due  to  the  wave  motion)  of  the  fluid  particles  is  sufficiently 
small  so  that  the  square  of  this  velocity  can  be  neglected  In  Bernoulli's 
equation  (Reference  1*.  page  1}  the  velocity  potential  of  fluid  motion  (for 
an  incompressible,  non-viscous  fluid)  due  to  a  source  (Reference  1,  page  404) 
located  below  the  free  surface  of  a  uniform  stream  of  infinite  depth  is  (Ref¬ 
erence  2,  page  3) 


•Feferecce*  are  listed  an  page  17 . 
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where  x,y,z  are  rectangular  coordinates,  z  positive  upwards  (the  undis- 
turbed  free  surface  is  the  xy-plane  of  Figure  1), 

f  is  the  depth  of  the  source  below  the  undisturbed  free 
surface, 

r2  is  equal  to  x2  +  y2  +  (z  +  f)2 

r2  is  equal  to  x2  +  y2  +  (2  -  f)2 

*(x,y,z)  is  the  velocity  potential  v  =  --%£>  v  =  - 

u.v.w  are  velocities  in  positive  x,  y,  z  directions, 

-c  is  the  uniform  stream  velocity, 

k0  is  equal  to  g/c2 

g  is  the  acceleration  of  gravity, 

2  is  the  strength  of  the  source  (a  source  of  strength  m  emits 
a  volume  tnm  per  unit  time) ,  and 

?  placed  before  an  integral  sign  means  that  the  Cauchy  prin¬ 
cipal  value  is  to  be  taken  (Reference  3,  page  '28). 


In  the  notation  used  the  source  is  con- 
*  aider ed  to  be  held  stationary  In  a  uni- 

/  form  stream.  This  is,  of  course,  equiv- 

/  a lent  to  having  the  source  move  with  a 

•0,0, "-^2  / - -a  uniform  velocity  through  a  stationary 

/  J  fluid.  It  is  of  interest  to  note  that 

/  V>.  »  the  first  two  terms  of  Equation  (1] 

L _ '-r- _ — — _ *  give  the  velocity  pot  ial  of  a  source 

in  an  unbounded  uniform  stream. 

If  a  source  and  equal  sink 

✓ 

_ _ Locmioo  or  so«c»  are  placed  a  given  distance  apart  on  a 

Jo,  a  -o  •— 

line  parallel  to  a  uniform  stream  in  an 
infinite  fluid,  the  resulting  fluid 
Figure  1  motion  is  that  for  the  flow  of  a  uni¬ 

form  stream  about  an  oval  shaped  body 

called  the  Rankine  Ovoid  (Reference  1,  page  4n).  If  the  source  and  sink  are 
near  the  free  surface  of  a  uniform  stream  the  motion  is  that  for  the  flow  of 
a  uniform  stream  (with  a  free  surface)  about  a  "distorted"  Rankine  Ovoid.  The 
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distortion  Is  small  when  the  source  and  sink  are  sufficiently  far  from  the 
free  surface  of  an  Infinitely  deep  stress.  Using  Equation  [1  ]  the  velocity 
potential  of  the  fluid  motion  about  the  "distorted"  Ranklne  Ovoid  with  source 
at  (a,o.-f)  and  sink  at  (-a,o,-f)  is 


#(x.y.z) 


cx  +  -5- 

ri 


+ 


a 

r4 


-  **’*0°  pf^sec2dd*  re~*</~,)coalk(x-a)cos»]cos(kysln»)  dk 

*  **  k-kosec20 

9 

-  4k0mjTsec2#e  sin[k0(x-a}sec*]cos{koyslndsec2#)d®  j2] 

+  .—■£??■  P  (Vc20d«  re'kl/t) cos(k(x+a}cos  eJcosQffisln  ^ 

9  J°  k-k0sec2e 

+  4k0m|sec2®e~*iU-*>**‘’z#  sin[k0(x+a)seca',cos{k<jysindsec2a}da 
Jo 

where,  in  addition  to  tens  defined  for  Equation  (l  ] 

r2  ■  (x-a)2  +  y2  +  (z+f)2 
r2  *  (x-a)2  +  y2  +  (z-f)2 
r2  *  (x+a)2  +  y2  +  (z+f)2 
rf  *  (x+a)2  +  y2  +  (z-f)2 

■file  first  three  terms  of  Equation  [2]  give  the  velocity  potential  of  a  3ource 
and  equal  sink  In  an  unbounded  uniform  stream. 

La gaily1 3  theorem4’5  nay  be  applied  to  obtain  the  moment  acting  on 
the  above  "distorted"  Ranklne  Ovoid.  This  theorem  states  that  the  forces  act¬ 
ing  on  a  body  whose  surface  Is  a  closed  stream  surface  of  the  fluid  motion, 
are  given  by  the  vectors 


P  =  -4jrpnq 


where  a  Is  the  strength  of  a  source  Internal  to  the  stream  surface, 

q  Is  the  resultant  fluid  velocity  vector  at  the  location  of  the 
source  due  to  all  other  sources  (Its  components  In  the 
x,  y,  z  directions  are  q^,  q,  q2),  and 

p  is  the  mass  density  of  the  fluid. 


Thus,  the  force  has  the  direction  of  -q  and  Its  line  of  action  passes  through 

the  point  at  which  the  3ource  is  located.  Each  tsrs  of  Equations  Ill  and  [2] 
say  be  considered  as  the  velocity  potential  due  to  a  certain  source  or  distri¬ 
bution  of  sources.  Applying  Lagally's  theoren,  the  vertical  force  acting  on 
the  "distorted"  Hankine  Ovoid  through  the  point  (a,o,-f)  is 

Z(a,o,-f)  «  -4irp«qz(a,o,-f )  [ 3] 

and  the  vertical  force  acting  on  the  body  through  the  point  (-a,o,-f)  is 

Z(-a.o.-f)  «  -4»p(Hi)q2(-a,o,-fi  [4] 

Since  there  are  no  other  internal  sources  for  this  body.  Equations  [3]  and  [4] 
give  the  only  vertical  forces  acting  on  the  body. 

Carrying  out  the  partial  differentiations  indicated  In  Equations  [3] 
and  (4],  the  following  expressions  are  obtained 
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Hence,  the  moment  acting  on  the  body  about  the  center  or  buoyancy, 

Kx  =  a'2'a.o.-f)  -  ia 

Hx  -  32*pa!n2k|J^sec4#e  tfkot0C  9sln[2sk03ec®]dfi  [5] 

(In  the  above  notatim  a  positive  moment  acts  to  raise  the  nose  of  the  body.) 
Making  the  substitution  t  =  tan#,  this  equation  nay  be  written  as 

Mx  *  32*j9aa?k|J" (l+t2)e~v*0<lkl*)  sin(2ak0Vi+t2 ]dt 

Por  Ranklne  Ovoids  with  sufficiently  large  length-diameter  ratios 


where  b  is  the  maximum  radius  of  the  Rankins  Cvoid.  With  this  substitution 
the  expression  for  Mx  becomes 

M  =  »pbVk0(2ako)e-V*o/  •~*/*0'*  O+t* )sin[2ak#Kl+t2  ]dt  [5*  ] 

*  *0 

THE  SECOND  APPROXIMATION 

Equation  [5]  gives  the  moment  acting  on  a  "distorted3  Rankin-;  Ovoid. 
To  obtain  a  closer  approximation  to  the  moment  acting  on  an  undistorted 
Rankine  Ovoid  placed  below  the  surface  of  a  uniform  stream  the  simple  distri¬ 
bution  of  a  source  and  equal  sink  on  the  axis  of  the  closed  stream  surface  rep¬ 
resenting  the  Rankine  Cvoid  must  Jbe  codified.  This  might  be  done  by  an  exten¬ 
sion  of  the  method  of  images  as  applied  by  Havelock  in  the  case  of  a  circular 
cylinder.  This  method  would  require  finding  the  image  system  of  sources  with¬ 
in  the  closed  stream  surface  due  to  the  source  system  above  the  free  surface. 
However,  instead  of  attempting  to  obtain  this  image  system  exactly  an  approx¬ 
imate  image  system  will  be  sought. 

It  has  been  shown  by  von  Kfoc/n*  that  for  a  body  of  revolution  with 
its  axis  parallel  to  a  uniform  stream  the  effect  of  superimposing  a  flow  per¬ 
pendicular  to  the  axis  may  be  obtained  approximately  by  a  suitable  distribu¬ 
tion  of  doublets  (Reference  6,  page  12)  along  the  axis  if  the  body  between  the 
limits  of  the  source-sink  distribution  which  defines  the  body  in  the  uniform 
stream.  The  doublets  are  oriented  so  that  their  axes  are  opposite  in  direc¬ 
tion  to  the  transverse  flow  and  their  strength  per  unit  distance  along  the 
axis  of  the  body  is 


ijh 
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a  «  -|r*w  161 

where  r  Is  the  radius  of  the  body  at  the  position  of  the  doublet,  and  w  is  the 
superimposed  transverse  velocity. 

Therefore  in  the  case  of  the  Rankine  Ovoid  aovlmg  below  a  free  sur¬ 
face,  the  effect  of  the  velocity  Induced  by  the  free  surface  way  be  accounted 
for  by  a  suitable  distribution  of  doublets  along  the  axis  of  the  s trees  sur¬ 
face  between  the  source  and  sink.  It  will  be  apparent  that  for  the  calcu¬ 
lation  of  aoaents  (see  below)  only  the  vertical  coaponent  of  the  Induced  veloc¬ 
ity  need  bo  considered.  This  vertical  velocity  is  obtained  fron  Equation  [2] 
by  evaluating at  points  along  the  line  between  the  source  and  sink.  Since 
the  vertical  velocity  does  not  change  very  rapidly  with  depth,  this  calcula¬ 
tion  of  the  vertical  velocity  is  probably  satisfactory.  For  Rankine  Cvoids 
with  fairly  large  length -diaweter  ratios  the  dlaseter  of  the  body  is  nearly 
constant  in  the  region  between  the  source  and  sink.  Thus,  as  a  further  simpli¬ 
fication,  the  radius  r  in  Equation  [6]  will  be  considered  constant  and  equal 
to  the  aaxlaua  radius  of  the  body.  The  desired  doublet  distribution  is  now 
given  by 

*-  -2**w  17] 

where  b  is  the  aaxlaiaa  radius  of  the  undlstarted  Rankine  ovoid,  and  w  is  the 
vertical  component  of  the  Induced  velocity  calculated  froa  Equati  on  [2]. 

An  application  of  legally' s  theorea  shows  that  a  body  whose  surface 
is  a  closed  streaa  surface  of  the  fluid  notion  experiences  a  mommt  if  the 
axis  of  any  Internal  doublet  is  noraal  to  the  direction  of  a  superimposed  uni- 
fora  streaa.  (The  doublet,  however,  causes  no  resultant  force  cm  the  body.) 
This  aoxent  is  given  by  (Reference  5,  page  13) 


N  «  -Kme*c 


18] 


where  p  is  the  aass  density  of  the  fluid, 
a  is  the  doublet  strength,  and 
c  is.  the  unifora  streaa  velocity. 

Hence,  froa  Equations  (7]  and  [8]  the  sonant  per  unit  length  for  the  Rankine 
Ovoid  (due  to  the  unifora  streaa  c)  has  the  algebraic  sign  of  the  vertical 
velocity  and  is  given  by 


H  »  2*pcb2w 


19J 
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and  the  total  moment*  on  the  body  Is 

Mg  *  2xpcb2£  w  dx 

where,  using  Equation  [2],  the  velocity  w  Is 

w  =  w(x.o.-f)  = 


no; 


X  •  X 

J-o 
x  «  -f 


=  2fn[(x-a)2  +  (2f)2}'  -2fm[(x+a)2  +  (2f)2] 


rV» 


4k0m  rf 


pf  seca<der^-^SijEi3L-?)co.s»3  ^ 

40  Ja  k-k„sec^9 


4k0n 


P^o^^or^^coslktx-fajcos#] 

k-k0sec2S 


ok 


dk 


in] 


+  Jrk^nnsec4fe"*/‘o',c*,sin[k<)(x-a)secS]d® 

— o 

-  4k^n0sec4®er*/ko**fZ#  sin[k0{x+a)  sec0]d0 

Tne  first  four  terms  of  Equation  [11]  give  that  part  of  w  which  Is  skew  sym¬ 
metric  about  x  =  0,  and  the  last  two  the  symmetric  part.  Since  .'he  integrals 
of  the  skew  symmetric  terms  vanish.  Equation  [10]  now  becomes 


Mg  -  i6jr/>cD2kj;ia  j^jJsec40e'J/to’‘cJ*sIn[kfl(x-a}sec$]d0  dx 


-jJJsec4®e'*/tc*,c2#sln(k0(x+a}secS]d®  dx 


[12} 


The  double  integrals  in  Equation  [12}  may  be  expressed  in  a  core 
convenient  fora.  Consider,  for  example,  the  second  of  these  integrals  and  let 


*7bc  lcsagltadlsf  1  Telocity  dae  to  the  vxre  syrtca  1»  sas.ll  coapsred  vlth  the  miforx  stress  Telocity 
c  Is  not  cccsldertd  la  eelcalatl^  Jfe.  la  addition  tie  asztcsl  actlc as  of  tie  sources  «-»*  doablets 
wlttla  the  closed  stress  ssrfWee  representing  the  body  gtw*  ao  ressltant  force  or  (Itftnoce  5, 

pe^es  5  to  5). 
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M  =  Kj  +  K2 

=  jrptVk0(2ak0)e'*/*°JV*/**‘l  (l+t2)sinl2ak#>'l+t2  ]dt 

+  4jrpb<c2k,e'*/*#£e"*/‘#‘Vl+t2(cosl2ak,V'l+t2]-l  )dt 
or 

K  =  npb*  c2k0e~ifk»[  e~tfkt*2  W+t*  |2akc^l+t2  sin[2ak,V:l+t2] 
4  4  cos  l2aka Vl^-t2  dt 

A  scinent  coefficient  is  given  by 

C  -  *  - 

M  ^pc2jrb2i 


where  p  is  the  nass  density  of  the  fluid, 
c  is  the  unifora  streaa  velocity, 

»rb2  is  the  naxiaua  cross-sectional  area,  and 
£  is  the  length  of  the  body. 

Then,  for  Rc.nkine  Ovolds  with  sufficiently  large  length -diameter 
noaent  coefficient  can  be  written  as 


e#f  f  e^1*  tf+t2!—  Wt2  sin  [j-ViTt2] 

+  4cos  [—  VUt2 

20s  f2*  l?2  I?2  J  l?2 


[14] 


[14>] 


M51 


ratios*  the 


H6} 


fjdt 


•For  the  present  the  terx  *saTflcl«ctly  large"  ma  10  or  greater.  Saskise  Oroide  with  other  lergti- 
dleaeter  ratio*  will  be  lsTeatlgated  is  a  later  report. 
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a  _  length  of  body  _  _5 
wnere  0  -  dlame»  er  c.f-y5Sy  ~  d  ’ 

depth  of  centerline  of  body  _  f_ 
a  *  “Jiaseter  of  body  =  d  ' 


distance  betveen  source  and  sink  2a  _ . 

«  - - riHirH - .  T  ,  and 


Pron  Equation  { 1 6  J  It  can  be  seen  that  approaches  zero  as  P  be¬ 
comes  either  very  small  or  very  .large .  The  two  cases  of  low  and  high  speeds 
nay  also  be  approximated  by  considering  the  free  surface  as  a  rigid  wall  for 
very  low  speeds  and  fcy  neglecting  gravity  forces  for  very  high  speeds.  In 
each  case  the  moment  coefficient  Is  zero. 

Writing  Equation  [l1* *  j  in  the  fora 


M  - 


d4 


£ ESI 


t/3 

C2 


IS 


sinf fJi  Vut2l 
-  °  l  C2  L  c2  J 


+  4cos 


dt 


it  is  evident  that  for  a  given  length,  deptn,  and  speed  the  moment  H  is  approx¬ 
imately  proportional  to  the  fourth  power  of  the  diameter  for  Ranklne  Ovoids 
with  length-diameter  ratios  of  10  or  greater.  The  approximation  depends  on 
(.  For  a  length-dlaaeter  ratio  of  10.5  «  is  0.95  and  as  the  length-diameter 
ratio  is  Increased  €  approaches  1  as  a  limit. 

Tne  integrals  appearing  above  nay  be  evaluated  either  by  numerical 
quadrature  or  by  expanding  the  integrand  in  an  infinite  series  md  integrating 
term  by  term.  The  series  expansion  is  rapidly  convergent  for  P  >  1  {I.e., 
for  large  speeds  of  advance)  while  for  P  <  1  It  is  more  convenient  to  use  a 
Gauss  type  quadrature  formula.  The  two  methods  are  described  in  the  Appendix. 

Pigures  5  and  4  indicate  the  values  of  C„  and  C„  «  C,.  +  C„  given 

”1  «  nj  Hi 

by  Equation  [16]  for  a  Ranklne  Ovoid  with  a  length-diameter  ratio  of  10.5  at 
submergences  of  1-5  and  2.8  diameters.  The  dotted  line  in  Pigure  4  shows  the 
values  of  obtained  from  the  approximate  formula  for  given  by  Equation 
(?l ]  of  tne  Appendix. 
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Figure  4  -  Moment  Coefficient  +  for  Ranklne  Ovoid, 
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CONCLUSIONS 

The  difference  between  the  first  and  second  approximations  for  the 
moment  acting  on  the  Ranklne  Ovoid  moving  below  a  free  surface  indicates  that 
It  is  essential  to  modify  the  singularity  system  representing  a  body  In  an 
infinite  fluid  in  order  to  account  for  the  effect  of  the  waves  formed  by  the 
motion  of  the  body  if  the  correct  moment  acting  on  che  body  is  to  be  obtained. 
To  the  order  of  approximation  considered  the  modification  of  the  singularity 
system  is  such  that  It  gives  rise  only  to  a  couple. 
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APPENDIX 

EVALUATION  OP  INTEGRALS  APPEARING  IN  TEXT 

As  mentioned  in  the  report  the  integrals  appearing  in  the  expres¬ 
sions  for  moments  can  be  evaluated  either  by  numerical  quadrature  or  by  expand¬ 
ing  the  integrand  in  an  infinite  series  and  integrating  term  by  term.  An 
efficient  method  of  numerical  quadrature  for  integrals  of  the  type  required  is 
given  by  the  Gauss-Christoffel  formula 

f  e'**  G(x)dx  —  5^A^G{x, }  (1?] 

1=0 


where  the  x^s  are  roots  of  the  Heraite  polynomials  of  the  order  and  the 
A^s  are  called  Christoff  el  numbers.  The  theory  of  this  method  of  numerical 
quadrature  is  given  in  References  8  and  9.  and  tables  of  the  x^'s  and  A^'s 
for  various  values  of  n  are  given  in  References  10  and  11. 

Since  the  x^s  are  symmetrically  spaced  on  either  side  of  x  =  0  and 
the  Aj^s  are  symmetric  about  x  =  0,  the  formula  can  be  used  when  G{x)  is  sym¬ 
metric  for  the  evaluation  of  integrals  of  the  fora 

f  e~,2G(x)dx 
Jo 

In  this  case  only  the  positive  x^s  of  Equation  1173  are  taken. 

As  an  application  of  this  formula  consider  Equation  [16]  for  the 
moment  C^-  Making  the  transformation 


P 


t 


the  resulting  expression  for  is 

CM  =  sf  e'p*G(p)dp 

n  jo 


Pi  8] 


where  s  = 


2 P?eVw 


u 


2a 

fi** 
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G(p)  «  $2sin$  +  4«(cos$-l),  and 

As  a  particular  example  consider  the  Ranklne  Ovoid  described  on 
page  10  snd  let  the  Proude  nu*be~  be  P  »  0.7  and  the  depth  of  submergence  be 
1-5  diameters.  The  constant  s  is  equal  to  0.0034808,  and  applying  the  quadra¬ 
ture  formula  (witn  n  *  4}  the  result  is* 


CM  =  0. 0034308 [x0G(p0)  +  XlG{pi)  +  -  •  •  +  A4G(pj] 

=  0.0034808  [0.36012  G ( 0 )  +  0.43265  G(0. 72355) 

+  0.088475  G(l.4o8b)  +  0.0049436  G(2.2666)  +  0.000039607  G(3-1910)J 
«  0.0034808  [-2.5525  -  7-3905  -  3.5838  -  0.037864  -  0.00047394] 

=  -0.047217 

Repeating  this  process  for  n  *  7  the  result  is  =  -0.047207 

No  general  expression  for  the  error  in  this  method  of  quadrature  is 
available.  However,  for  the  case  considered  here,  a  check  can  be  obtained 
since  the  integral  can  be  expressed  as  an  alternating  series  from  which  the 
correct  answer  nay  be  determined  to  any  desired  accuracy.  This  series  expan¬ 
sion  is  discussed  in  the  next  few  paragraphs. 

The  corresponding  series  expansion  for  is  obtained  by  expanding 
tne  integrand  in  an  infinite  series  and  then  (since  the  series  is  uniformly 
convergent)  inverting  the  order  of  integration  and  summation.  Thus,  substi¬ 
tuting  the  series 


n+r 


Kl+t*  sinj^  yi+tzj 


■  2i-f 

n=l 


(2n)i 


0'  A  ari 


t. 


-ror  FeJ>rtace  1C. 


15 


In  Equation  [16],  the  result  Is 


*n_ .» 


20s  P2 


Havelock  ha3  observed  (Reference  12,  page  284}  that  integrals  of  the  form 


L2».,  =j0e"’’  (1+t;2)  dt 


where  n  is  an  integer ,  can  be  expressed  in  terms  o  modified  Bessel  functions 
of  tne  second  kind  of  orders  zero  and  one.  One  procedure  for  doing  this  de¬ 
pends  on  the  formula* 

wfe-W^-dt  =  fe’^O+t2)12*1"  dt 

Jo  *•  Jo 


2^joe-'2(i+t2p^  dt 


For  2n+i  =  3  this  formula  gives 


L  =  f*e‘,,'2(l+t2)s/2  dt  =  f  e*n,1(l+t2)I/2dt  -  f  e'^Cl+t2)' ‘^dt 
3  Jo  r>  Jo  2*1  Jo 


1  +  1?  T  .  1_T 

TJ  \  27}  L-1 


ana  continuing  the  procedure. 


=  2+2  r  3  L 

7]  "3  271  1 


3  7} 


l  _  3+2.  l  -  i 

"Z«*l  ~  jj  27}  2«-3 


•T:is  rc-octicr  ror=uli  vas  scour.  to  tre  author  by  Dr.  J.V.  Wrench,  Jr.  It  aav  he  derived  by  tat! 
t: *■  cr'.v--.  ive  v:tr.  respect  to  t  of  Z»+l 

t(l+t2 T~  e'15'2 

rearrir^ir terns.  and  t:  er.  Integrating  hetveen  the  Units  0  and  ». 
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Hence,  from  Lx  and  L.j  all  the  other  L's  nay  be  determined.  The  transforma¬ 
tion  t  -  sinh/i  applied  to  thb  integrals  for  Lx  and  L.x  transforms  them  into 
known  expression*  for  modified  Bessel  functions  of  the  second  kind  K(  and  1L 
(Reference  13,  page  181).  The  result  is 


« 

v  >“[*.(!)  * «.(?)] 


i  i 


Since  the  functions  X#(x)  and  Kx(x)  have  been  tabulated. >*’H  the  L- 
f unctions  can  be  determined.  The  expression  for  Cg  now  becomes 


n-i 


120] 


For  the  same  example  considered  on  pages  13  and  14  the  sum  of  the 
first  11  terms  of  this  series  is  -0.0472070  and  the  12th  term  is  -0.0000007- 
Hence  to  4  significant  figures  the  value  of  Cg  is  -0.04721.  This  shows  that 
the  quadrature  formula  with  n  *  7  gives  the  correct  answer  to  4  significant 
figures  and  for  n  -  4  gives  an  answer  that  is  only  in  error  by  1  in  the  fourth 
significant  figure.  Since  3  significant  figures  give  sufficient  accuracy  for 
any  application,  the  quadrature  formula  with  n  «  4  was  used  in  calculating  the 
results  given  in  Figures  3  and  4. 

Noting  that  the  term  for  n  *  2  is  zero  in  the  series  [20]  it  appears 
that,  for  sufficiently  large  F,  the  first  tern  of  the  series  may  be  a  good 
estimate  of  Cg.  Figure  4  compares  the  correct  Cg  with  the  values  obtained 
from  the  first  term  of  Equation  [20]  for  a  Rankin e  Ovoid  with  a  length-diameter 
ratio  of  10-5  submerged  1.5  diameters  and  2.8  diameters.  The  figure  shows 
that  for  Froude  numbers  F  greater  than  1  the  first  term  of  Equation  [20]  is  a 
good  approximation  to  Cg. 

Thus  for  values  of  the  parameters  not  too  different  from  those  in 
the  example,  the  approximate  value  of  the  moment  coefficient  for  Rankine 
Ovolds  is 

c« = ’  5?  "***  [wK*@) +  (i+w)ki(?)]  1 


where 


u  = 


?q 
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